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L'invention concerne une methode pour discriminer des signaux microsismiques 
induits par l'exploitation d'un gisement souterrain, parmi des signaux sismiques emis dans 
le cadre d'operations de surveillance sismique 

Pour suivre revolution d'un gisement souterrain en cours d'exploitation, il est connu 
d'implanter a poste fixe des sources sismiques et des recepteurs sismiques, en surface ou 
dans des puits, qu'il s'agisse de puits d'exploration ou d'exploitation du gisement et de 
proceder a intervalles de temps fixes a des operations d'enregistrement des signaux 
sismiques emanant de la formation. II peut s'agir de signaux sismiques induits par l'activite 
d'exploitation du site ou de signaux sismiques renvoyes par les discontinues geologiques 
15 de la formation, en reponse a remission d'ebranlements dans le sol au moyen d'une ou 
plusieurs sources sismiques. La comparison des enregistrements effectu6s a des instants 
differents de l'exploitation du site, fournit des donnees sur son Evolution. 

Differents systemes de surveillance sismique a long terme sont decrits par exemple 
dans les brevets EP 591 037 (US 5 461 594), FR 2 593 292 (US 4 775 009), FR 2 728 973 
20 (US 5 724 311) ouFR 2 775 349. 

Par le brevet EP 748.457 (US 5 724 311) de deux des demandeurs, on connait un 
systeme de surveillance permanente permettant d'obtenir une parfaite reproductibilite des 
conditions operatoires dans des operations repetitives de surveillance sismique d'une zone 
souterraine traversee par au moins un puits ou forage, et notamment d'un reservoir 
25 souterrain de stockage de gaz. Ce systeme comporte, installs a poste fixe, un ou plusieurs 
ensemble(s) de recepteurs sismiques (enterr6s en surface ou places dans un ou plusieurs 
puits), plusieurs sources sismiques repetitives (enterrees ou en surface) ainsi que d'un 
reseau de liaison permanente pour l'alimentation selective de ces sources en energie. Une 
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station centrale t€lecommande s61ectivement chacune de ces sources, et enregistre les 
signaux sismiques emanant de la zone souterraine en reponse aux ondes sisrniques 
transmises s61ectivement dans le sol par les sources. 

Grace k cet ensemble de sources & poste fixe dont le couplage avec les terrains 
5 environnant reste stable, et & ce r€seau d f alimentation au moins en partie enterre et dont 
l'aire d'emprise en surface est r6duite, on peut mener toute une serie d'operations sismiques 
de surveillance de longue dur6e dans des conditions op6ratoires stables, sans risque 
d'incompatibilit6 avec les activites du chantier d'exploitation. 

Par le brevet FR 2.775.349 (US 6 182 012) des deux mdmes demandeurs, on 
10 connait une methode et un dispositif de surveillance permanente d f une formation 
souterraine dans laquelle on forme des cavit6s pour des sources sismiques (de preference 
assez prof ondes pour qu'elles soient couplees avec la formation au-dessous de la' .zone 
alt6r£e (WZ)), on fore au moins un puits soit au-dessous de chacune des sources ou au 
voisinage imm6diat de chacune d'elles. Dans chacun des puits ainsi for6s, on descend une 
15 antenne constitute d'une plurality de r6cepteurs sismiques (g6ophones, hydrophone, gtc.) 
que Ton couple avec la formation environnante. Avec ce dispositif, on realise des cycles 
d'emission-r6ception avec emission Sondes sismiques dans la formation par au moins fine 
source et acquisition des signaux capt6s par les r6cepteurs sismiques, en reponse aux ondes 
renvoyees par la formation. Les r6cepteurs peuvent par exemple Stre associes avec un tube 
20 de transfert de fluides descendu dans le puits, dans le but de relier une zone souterraine h 
un appareil d'exploitation en surface et couples avec la formation environnant le puits. 

Les differentes sources du dispositif de surveillance sismique peuvent etre 
actionnees successivement, en m6nageant entre les declenchements un intervalle de temps 
suffisant pour recevoir les ondes renvoy6es par la zone investigu6e. On peut aussi utiliser 
25 plusieurs sources sismiques 6mettant les memes signaux que Ton declenche simultanement 
pour augmenter la puissance 6mise. 

Par le brevet FR 2.805.051 des demandeurs, on connaft une m6thode et un syst&me 
de surveillance sismique d'une formation souterraine. Comme schematises aux Fig.l &4, le 
syst&me comporte par exemple un reseau d'antennes sismiques 2 constitutes chacune d'un 
30 ensemble de capteurs sismiques 4 disposts^ intervalles rtguliers le long d'un puits 3 for6 
dans le sol. Ce reseau peut etre r6gulier comme sch6matis6 sur la Fig.2, ou irr6gulier. Les 
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capteurs sont par exemple des gdophones mono-directionnels orients verticalement ou 
multi-axes (triphones) et/ou des hydrophones. A proximite de chaque antenne 2, est 
dispos6e une source sismique 5. Comme sources, on utilise avantageusement des vibrateurs 
de type pi6zo61ectrique, tels que d6crit dans le brevet FR 2.79 1.780 (US 6.338.394) aux 

5 noms conjoints des demandeurs, qui sont installs a demeure au voisinage immediat de 
chaque antenne 2. A chaque antenne peut etre associe un module 61ectronique de controle 
et d'enregistrement 6 (Fig.l). Ces modules 61ectroniques 6 peuvent etre reli6s & un 
laboratoire central de commande et de synchronisation 8. Suivant un autre mode de mise en 
ceuvre, schematise sur la figure 4, les antennes 2 sont reli6es directement & un laboratoire 

10 central 8 regroupant tous les modules 61ectroniques individuels de controle et 
d'enregistrement 6 (Fig.l). 

Des ondes sismiques sont 6mises dans la formation souterraine par une ou plusieurs 
sources sismiques (vibrateurs) 5. Dans le cas de sources multiples, celles-ci 6mettent 
simultanement ; elles sont pilot6es par des signaux orthogonaux les uns relativement aux 

15 autres, de maniere & former un signal vibratoire composite. Les vibrateurs fonctionnent en 
permanence ou tout du moins sur une dur6e assez longue pour obtenir un rapport signal au 
bruit suffisant. Les ondes sismiques gen6r6es par les sources sismiques 5 se propagent vers 
le bas (ondes descendantes 9). Ces ondes incidentes sont tout d'abord enregistr6es par les 
recepteurs 4 de chaque puits 3. Les ondes renvoy6es par les discontinues de la zone 

20 (interfaces sismiques) se propagent vers le haut. Ces ondes montantes 10 sont enregistrees 
aussi par les differents recepteurs 4. De la sorte, les ondes montantes et descendantes sont 
superpos6es sur les sismogrammes. On applique aux enregistrements les traitements usuels 
bien connus des gens de Tart convenant & ceux obtenus par la technique dite de PSV 
(Profils Sismiques Verticaux). Ensuite, on discrimine les contributions respectives des 

25 sources sismiques au signal vibratoire composite et on reconstruit les sismogrammes 
Equivalents k ceux que Ton obtiendrait en actionnant separ6ment les sources sismiques. 
Comme signaux orthogonaux, on utilise par exemple des signaux sinusoi'daux de 
frequences differentes les unes des autres aussi bien par leurs composantes fondamentales 
que par leurs harmoniques respectives, ou des signaux formes i base d'ondelettes, de 

30 polynSmes de Legendre ou de s6ries aldatoires, etc. 

Dans le cas notamment ou les signaux orthogonaux 6mis sont des sinusoides, on 
realise par exemple la discrimination des contributions respectives des sources sismiques, 




par determination de P amplitude et de la phase du signal vibratoire composite aux 
frequences fondamentales des signaux de pilotage appliques aux sources sismiques. Cette 
discrimination comporte par exemple une ponderation des signaux enregistres par un 
facteur de ponderation (ou d'apodisation) en cloche et une determination de P amplitude et 
5 de la phase du signal composite, ou bien encore une selection par transform6e de Fourier, 
de raies du spectre complexe associees respectivement aux differents signaux ponder£s. La 
reconstruction des sismogrammes correspondant specifiquement aux diff6rentes sources 
sismiques est effectu6e par exemple en appliquant, apres leur separation, une transformee 
de Fourier inverse, aux raies associees respectivement aux differents signaux ponderes. 

10 Par ailleurs, il est connu que P exploitation des reservoirs induit des variations de 

contraintes mecaniques dans le reservoir et dans sa couverture. Ces variations de contrainte 
ont pour effet de generer des signaux sismiques se propageant dans les terrains. Les 
recepteurs des differentes antennes 2 vont enregistrer ces signaux (ondes montantes 11). 

Comme les sources sismiques 5 sont situees beaucoup plus pres des recepteurs 4 
15 des antennes 2 que du reservoir et que renergie qu'elles emettent est bien plus intense que 
renergie des signaux microsismiques induits genets dans le reservoir, il en resulte que les 
signaux microsismiques sont gen6ralement masques et indecelables. 

La methode selon Pinvention 

La methode selon invention permet de separer des signaux de microsismicite 
20 induite et des signaux sismiques multi-sources dans le cadre d f operations de la surveillance 
d'une zone souterraine en cours Sexploitation > 

Elle comporte Pemission d' ondes sismiques dans la formation en couplant avec elle 
une ou plusieurs sources sismiques. Dans ce demier cas, celles-ci emettent simultanement 
des signaux orthogonaux les uns relativement aux autres, de manifere k former un signal 
25 vibratoire composite, la reception des signaux renvoyes par la formation en reponse a 
pemission des ondes sismiques, Penregistrement des signaux regus par au moins un 
capteur sismique et le traitement des signaux enregistres pour separer les contributions 
respectives des sources sismiques aux signaux regus et reconstruire les sismogrammes 
equivalents a ceux que Pon obtiendrait en actionnant separement les sources sismiques. 
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On separe les signaux de microsismicite induite (sismique passive) des signaux 
sismiques 6mis essentiellement en isolant la contribution de chacune des sources par 
comparaison avec un modele spectral de reference aux frequences fondamentales 6mises et 
& leurs harmoniques respectives, et en reconstituant dans le domaine des temps les signaux 
5 de microsismicite induite et les signaux sismiques 6mis. 

La contribution spectrale des signaux de microsismicite au spectre des signaux 
re9us est obtenue par exemple en retranchant les valeurs Samplitude et de phase associ6es 
au module spectral de reference, de 1' amplitude et de la phase du spectre associe aux 
enregistrements. 

10 Suivant un mode de mise en oeuvre, comme modele spectral de reference, on utilise 

par exemple un module courant forme par mise k jour d'un modele spectral prenant en 
compte la contribution de cycles d'enregistrement precedents. 

Le module spectral courant peut etre forme, en determinant une valeur moyenne des 
spectres de frequence formes a partir d'enregistrements (anterieurs et/ou posterieurs) 
15 obtenus pour la meme source et les mSmes frequences. 

Suivant un autre mode de mise en oeuvre, on forme le modele spectral courant, en 
determinant une valeur m£diane des spectres de frequence formes & partir d'enregistrements 
anterieurs obtenus pour la mSme source et les mSmes frequences. 

Suivant un autre mode de mise en oeuvre, on forme le modele spectral courant, par 
20 extrapolation ou interpolation k partir du spectre de frequence & partir de valeurs spectrales 
voisines. 

Suivant un autre mode de mise en oeuvre, on forme le modele spectral courant en 
utilisant la raie obtenue k cette frequence. 

Suivant un mode de mise en oeuvre, la discrimination de signaux de microsismicite 
25 induite parmi des signaux emis dans le cadre d'operations de surveillance sismique active 
est obtenue par les etapes suivantes : 

a) pour chaque enregistrement n d'un cycle p d'enregistrement, on calcule les 
contributions respectives des differentes sources aux frequences fondamentales ; 




b) on calcule ensuite le rapport de la contribution k un modele spectral courant 
forme par actualisation d'un modele spectral precedent & partir de frequences emises durant 
l'enregistrement precedent ainsi que de leurs harmoniques ; 

c) on d£duit la partie de l'enregistrement n du cycle p, que Ton peut associer aux 
5 operations de surveillance sismique active ; 

d) on deduit la partie de l'enregistrement n du cycle p, relevant proprement de 
l'activite microsismique passive ; et 

e) on forme par inversion dans le domaine des temps, des contributions spectrales 
respectives de chacun des sources sismiques aux frequences fondamentales et k leurs 

10 harmoniques, apres completion d'un cycle de mesures, les sismogrammes que Ton peut 
associer aux operations de surveillance sismique active ; et 

f) on forme par inversion dans le domaine des temps, k partir de la partie relevant 
de l'activite microsismique passive, les signaux microsismiques sous-jacents contenus dans 
les enregistrements. 

15 La contribution spectrale est obtenue en multipliant une fonction de transfert entre 

une ondelette caracteristique de la source et un sismogramme associe a un recepteur donne, 
par une ondelette caracteristique de la source. Cette fonction de transfert peut etre 
continuellement mise k jour au cours d'un cycle courant k partir d'une estimation faite au 
cours d'un cycle precedent et d ! une estimation initiale faite au cours du cycle courant par la 

20 relation. 

Presentation sommaire des figures 

Les caracteristiques et avantages de la methode selon Tinvention, apparaitront plus 
clairement a la lecture de la description ci-apres d f un exemple non limitatif de mise en 
oeuvre, en se referant aux dessins annexes oh : 

25 - la Fig.l montre schematiquement le dispositif de surveillance mis en place sur le terrain 
pour surveiller sismiquement une formation souterraine ; 



la Fig.2 montre un mode de disposition des unites d'emission-reception sur le terrain ; 




- la Fig.3 montre schernatiquement une de ces unites d'emission-reception coraprenant 
une source sismique et une antenne de reception ; 

- la Fig.4 montre une variante du dispositif de surveillance de la Fig.l ; et 

- la Fig.5 montre des exemples de signaux et spectres de frequence mentionnes dans le 
5 cours de la description. 

Description detailiee 

La methode permet done de r6aliser des operations de surveillance sismique d'une 
zone souterraine en utilisant un ensemble de capteurs sismiques et une ou plusieurs sources 
sismiques vibratoires. Dans le cas d'utilisation de plusieurs sources simultanSes, celles-ci 
10 sont actionn6es simultandment par des signaux k des frequences differentes choisies de 
fagon que Ton puisse discriminer les contributions de chaque source sur les sismogrammes 
obtenus & partir des signaux regus et enregistres. 

Ceci est realise d'une fagon g^nerale par le pilotage des diff6rentes sources par des 
sinusoides de frequences differentes et par utilisation de techniques de calcul numerique 
15 connues telles que la transformee de Fourier inverse pour separer les contributions aux 
sismogrammes obtenus des differentes sources sismiques. 

La distorsion ne peut etre negligee avec des sources m6caniques. En mSme temps 
qu*une frequence fi, une meme source Si 6met les frequences 2f u 3£ ... nfj. En 
consequence, si fi et fj sont les frequences respectives des deux sources Si et Sj de 
20 F ensemble de sources, il convient si on veut separer leurs contributions respectives non* 
seulement que fi * fj mais aussi que f j ^ 2fj, f\ * 3fj , . . . fi ^ nfj . 

Comme deja decrit dans le brevet FR 2.805.051 precite, si on emet le signal 
composite P t constitu6 de la somme de N sinusoides {£, Ai, <S>, } avec 1 < i < N, toutes les 
frequences etant contenues dans une bande spectrale comprise entre deux frequences 
25 limites f b et f f , le sismogramme T t observe au point R aura pour transformee de Fourier h la 
frequence fj„ le nombre de module Aj et de phase 3> f egaux a 1' amplitude et 2t la phase de 
la sinusoide TV On peut ainsi, en emettant successivement toutes les sinusoides de 
frequences f b & f f , reconstruire par transformee de Fourier inverse le sismogramme Tt. 




Dans Ie cas ou, par exemple, toutes les amplitudes Ai sont egales a 1 et toutes les 
phases <E>, sont nulles, le signal P t obtenu est tres proche de celui resultant de 1' auto- 
correlation d'un signal a frequence glissante variant dans rintervalle [f b - f f> ] (sweep), 
utilise couramment en vibrosismique. D'apres la theorie de la transformed de Fourier 
5 discrete, bien connue des gens de l'art, si l'on desire ecouter la source Si pendant le temps 
t e , V increment de frequence entre les sinusoTdes est A/ = 1/te et le nombre de sinusoides 
necessaires est Nf = (f r - fb)te- 

On peut done exciter simultanement N vibrateurs SI, S2, Si,... S N installes sur le 
terrain au moyen de signaux vibratoires de frequences telles que chaque source est excitee 

10 successivement par chacune des N f sinusoides ci-dessus a chaque instant, sous reserve que 
les frequences respectives des sinusoides emises a un meme instant par les differentes 
sources sismiques sont toutes differentes les unes des autres. On r6alise ainsi avec chaque 
vibrateur N f enregistrements correspondant aux N f frequences selectionriees dans 
l'intervalle de balayage. La separation des signaux recus par les capteurs sur le terrain, en 

15 reponse a remission simultanee des differents signaux est ainsi obtenue par selection de la 
raie a la frequence appropriee. "i* ■ 

On repete ainsi des cycles de Nf enregistrements. L'enregistrement En d'un cycle 
d' acquisition donn6 contient un ensemble de frequences fondamentales et harmoniques qui 
ne depend que de n. Partant de cette observation, la methode que Ton va appliquer pour 
20 d&ecter les effets de la rnicrosismicite, repose sur le constat v6rifi6 dans la pratique, que les 
variations d' amplitude et de phase observers a la frequence f et a 1' instant t sont ties peu 
differentes de celles observ6es a la fr6quence / - A/ et a 1' instant t - At . Cela permet de 
former un modele spectral cumuiatif construit a partir d'enregistrements de cycles 
precedents aux mSmes frequences mais forme done sur un temps long pendant lequel les 

25 caract6ristiques des signaux sismiques emis sont susceptibles de varier, et adapts en 
permanence a partir d'enregistrements obtenus pendant une periode suffisamment courte 
pour que les variations de signal puissent etre n6gligees L'efficacit6 de la m6thode 
propos6e est optimale quand les frequences successives Smises par une meme source 
s'incrementent regulierement. On pourra, par exemple, alterner des cycles en frequences 

30 montantes et en frequences descendantes. 
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Pour parvenir dans la pratique k la mise en Evidence de cette contribution spectrale 
des signaux microsismiques sous-jacents, on va proceder comme suit 

Notations : 

D fi,n, va d6signer la frequence fondamentale 6mise par la source i pendant 
l'enregistrement n (cette frequence est independante du cycle) ; 

i 

Pour simplifier les notations des grandeurs suivantes, Tindice designant le num6ro 
du recepteur est omis. Quel que soit celui-ci, les grandeurs sont calcul6es de la meme 
manifere. 

a Rp, n va designer l'enregistrement de rang n du cycle p ; 

H Cp,i, n > la contribution de la source i k l'enregistrement n du cycle p ; 

a M PiI1> le module initial de la partie active de l'enregistrement n du cycle p. Ce 
modele ne contient que les frequences emises ainsi que leurs harmoniques pendant 
l'enregistrement n ; 

n Ep, n , le rapport entre les contributions actives de l'enregistrement n du cycle p et 
le modele initial de la partie active de l'enregistrement n du cycle p ; 

a A px , la partie active de l'enregistrement n du cycle p ; 

B Pp,n» la partie passive de l'enregistrement n du cycle p (microsismique) ; 

B k le coefficient d'actualisation des ecarts (typiquement de l'ordre de 15 k 

25%) ; 

D h le coefficient d'actualisation des modeles (typiquement de l'ordre de 5 k 10%). 

Compte tenu de ces notations, la m6thode consiste dans un premier temps, pour 
chaque enregistrement n du cycle p, au calcul des contributions respectives C Pt f tn des 
diff6rentes sources aux frequences fondamentales ainsi qu'k leurs harmoniques par la 
methode d6crite dans le brevet PR 2.805.051 pr6cite. Du fait de la pr6sence 6ventuelle 
d'6nergie microsismique, la contribution des difterentes sources est susceptible d'Stre 
bruit6e. 
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On peut ameliorer cette estimation . en associant les contributions de tous les 
rScepteurs d'un meme enregistrement C Pf i ia . Ces contributions peuvent etre consid£r6es 
comme le produit de deux facteurs dont Tun ne depend que de la source (et non du 
recepteur) et T autre du recepteur et non de la source. 

Si Tindice r represente le recepteur, 
P = S T . 

^p 9 t 9 n,T p,i,n A p 1 i,n,r 

T ne d6pend de p, i et n que par la variation des caract&istiques du r6cepteur et de 
son environnement avec le temps. Ces variations sont tres lentes et peuvent Stre negligees 
ou estim6es de fagon plus precise. 

S et T sont d£termin6s a un facteur pres. On peut ainsi choisir pour S P 4, n la moyenne 
des ondes descendantes. Cette moyenne est obtenue par un calcul bien coniiu de ceux qui 
traitent les PSV (mise & plat des ondes descendantes et moyenne). Les Tp^u/isont alors les 
fonctions de transfert entre Tonde descendante et le sismogramme. 

Dans le cas ou Ton desire prendre en compte la variation de T, on peut utiliser par 
exemple le coefficient d'actualisation h pour mettre h jour T P 4, n , r & partir *de T p> i (1 ,-i ir en 
procedant de la mani&re suivante : 

Si Co P( i,n t r d6signe l'estimation initiale de la contribution C et To p ,i >D| r, celle du 
facteur recepteur, on calcule S par mise k plat et moyenne des ondes descendantes puis 
Top,^ par la relation : 

To • = COpj -° 

°p,i,n 

On obtient l'estimation finale du facteur recepteur par : 
T p , u>r = (l-h)T Pf u_ u +hTo piin>r 

et la contribution finale de la source i au cycle p de T enregistrement n par : 

p — o T 
^p.i.n.r °p,i,n i pj,a,r 
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On calctde ensuite le rapport Epj, de 1' eriregistrement courant au module spectral 
par la relation : 

E =(1- k)E p + 

Ce calcul necessite la connaissance de M P|11 dont le calcul d6crit ci-dessous est 
5 effectue lors du cycle p-1. II comporte un lissage sur les frequences pr6cedentes de fagon & 
limiter l'effet de T6nergie microsismique dans son estimation. La formule proposee 
n'utilise que les enregistrements pr6c6dents et autorise done le temps r£el. Dans les cas 
difficiles, on peut remplacer cette somme par une moyenne ou une m^diane courante. H 
faut alors travailler en temps diff6r6 (te temps n&essaire & r£aliser la seconde moitie des 
10 enregistrements utilises dans le calcul de la moyenne). 

Si M p+M repr6sente le module initial de la partie active de l'enregistrement n du 
cycle p+1, alors : 

Ce modele consiste essentiellement en une moyenne pond6r€e de la contribution 
15 courante et des contributions pr6c6dentes. La formule utilis6e autorise la variation du 
module : la contribution du cycle p-q au modele utilise pour le cycle p est multipli6e par le 
coefficient (l-h) q qui tend vers z6ro d'autant plus vite que h est voisin de 1. Cela permet au 
module M pn de suivre les variations saisonnieres, et aux rapports E pn de ne refl6ter que les 
variations journalieres.c) Connaissant le modele spectral M Ptn ainsi que le rapport E P( n de la 
20 contribution au module precedemment etablie, on en deduit la partie A p , n de 
Tenregistrement n du cycle p, relevant de la sismique active Partie active : 

De la meme fagon la partie P p , n de l'enregistrement n du cycle p, relevant 
proprement de Tactivit6 microsismique (passive) se calcule par la relation : 
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Pp,D ~ Rp.n - A P|D 




Les contributions spectrales respectives C p ,i >n de chacune des sources Si aux 
frequences fondamentales et & leurs harmoniques, aprfes completion d'un cycle de mesures 
fournissent par FFT inverse, les sismogrammes (dans le domaine des temps) utilises dans 
le cadre de r exploitation active. 

5 Par FFT inverse, & partir de la partie P p , n relevant de contributions passives, on 

obtient 1' estimation de Factivite microsismique sous-jacente contenue dans 
renregistrement R p , n . 

La mise k jour du rapport Ep, n entre la contribution courante et le modele spectral 
Mp, n , est effectuee de fa?on relativement rapide de fa§on a permettre des variations rapides 
10 du signal. La mise & jour du modele initial s'effectue de fagon beaucoup plus lente ; en fait 
le modele initial pourrait ne pas etre mis k jour. 

L'interet de ces mises & jour est d'optimiser r estimation du signal actif.jquand ses 
variations ne peiivent Stre controiees (variations saisonnieres et journalieres). Dans les cas 
les plus favorables, elles ne sont pas necessaires et 1' estimation de la partie active se 
15 reduira au modele initial. 

Une autre simplification possible est de n^gliger le bruit sismique aux frequences 
emises (et, si besoin, aux harmoniques de ces frequences). On ne se pr£occupe plus 
d'estimer la partie active et on se contente, pour obtenir la partie passive, d'annuler ou 
d'interpoler les raies correspondant aux frequences emises (et, s'il le faut, celles de leurs 
20 harmoniques). 

Dans un premier temps, on supprime les raies harmoniques du spectre de frequence 
de chaque signal enregistre et on interpole 1' amplitude et la phase aux valeurs de raie h la 
frequence de la source ; et puis on reconstitue par inversion renregistrement non affecte 
par les emissions des sources contenant eventuellement un signal microsismique. Ce mode 
25 de mise en oeuvre convient plutdt dans le cas d'un nombre de sources peu eieve et/ou 
presentant une variabilite importante du signal avec le temps. 
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REVENDICATIONS 

1) Methode de surveillance sismique active d'une formation souterraine (1) 
permettant une separation de signaux de microsismicite induite de signaux sismiques emis 

5 dans le cadre d'operations de surveillance sismique active d ! une zone souterraine en cours 
Sexploitation, comportant la realisation de cycles d'enregistrement sismique avec Emission 
d'ondes sismiques dans la formation en couplant avec elle une ou plusieurs sources 
sismiques (5) qui, dans ce cas, 6mettent simultan6ment des signaux orthogonaux les uns 
relativement aux autres, de manidre & former un signal vibratoire composite, reception des 

10 signaux renvoyes par la formation en reponse a remission des ondes sismiques, 
enregistrement des signaux re<jus par au moins un capteur sismique (2) et traitement des 
signaux enregistr^s pour s6parer les contributions respectives des sources sismiques aux 
signaux re^us et reconstruire les sismogrammes equivalents h. ceux que Ton obtiendrait en 
actionnant separ^ment les sources sismiques, caracterisee en ce que Ton separe sur les 

15 enregistrements les signaux de microsismicite induite des signaux sismiques resultant des 
operations de surveillance active, en isolant leur contribution par comparaison avec un 
modele spectral de reference en tenant compte des contributions spectrales de chaque 
source (Si) aux frequences fondamentales emises et a leurs harmoniques respectives, et en 
reconstituant par inversion dans le domaine des temps les signaux de microsismicite. 

20 2) Methode selon la revendication 1, caracterisee en ce que la contribution spectrale 

des signaux de microsismicite au spectre des signaux regus est obtenue en retranchant les 
valeurs d'amplitude et de phase associ£es au modele spectral de reference, de 1' amplitude et 
de la phase du spectre associe aux enregistrements. 

3) Methode selon la revendication 1 ou 2, caract£ris6e en ce que le module spectral 
25 de reference est un module courant form6 par mise & jour d'un modele spectral precedent en 

tenant compte de la contribution spectrale de cycles d'enregistrement precedents. 

4) Methode selon Tune des revendications precedentes, caracteris6e en ce que Ton 
forme le modele spectral courant en determinant une valeur moyenne des spectres de 
frequence formes & partir d'enregistrements anterieurs et/ou posterieurs obtenus pour la 

30 meme source et les memes frequences. 
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5) Methode selon Tune des revendications pr6c€dentes, caract6risee en ce que Ton 
forme le modele spectral courant en determinant une valeur mediane des spectres de 
fr6quence formes k partir d'enregistrements anterieurs et/ou posterieurs obtenus pour la 
meme source et les memes frequences. 

5 6) Methode selon Tune des revendications pr6c6dentes, caract€ris6e en ce que Ton 

forme le modele spectral courant par extrapolation ou interpolation k partir du spectre de 
frequence & partir de valeurs spectrales voisines. 

7) Methode de surveillance sismique active d'une formation souterraine permettant 
une discrimination de signaux de microsismicite induite parmi des signaux 6mis dans le 

10 cadre depurations de surveillance sismique active d'une zone souterraine en cours 
Sexploitation, comportant la realisation de cycles d'enregistrement sismique avec emission 
d'ondes sismiques dans la formation en couplant avec elle N sources sismigues (Si) 
emettant simultanement et pilotes par des signaux orthogonaux les uns relativ^ment aux 
autres, de maniere & former un signal vibratoire composite, reception des signaux^renvoyes 

15 par la formation en reponse h V emission des ondes sismiques, enregistrement d%s signaux 
re?us par des moyens de reception sismique et traitement des signaux enregistres pour 
separer les contributions respectives des sources sismiques aux signaux .^regus et 
reconstruire les sismogrammes equivalents & ceux que Ton obtiendrait en actionnant 
separement les sources sismiques, caracterisee en ce quelle comporte les etapes suivantes : 

20 a) pour chaque enregistrement n d ! un cycle p d'enregistrement, on calcule les 

contributions respectives (C Pii( „) des differentes sources aux frequences fondamentales ; 

b) on calcule ensuite le rapport (Ep, n ) de la contribution a un modele spectral (M P( n) 
courant forme par actualisation d'un modele spectral precedent (M p , n .i) a partir de 
frequences emises durant 1 'enregistrement (n-1) precedent ainsi que de leurs harmoniques ; 

25 c) on deduit la partie (A p ,„) de T enregistrement n du cycle p, que l'on peut associer 

aux operations de surveillance sismique active ; 

d) on deduit la partie (P p>n ) de 1' enregistrement n du cycle p, relevant proprement de 
l'activite microsismique passive ; et 
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e) on forme par inversion dans le domaine des temps, des contributions spectrales 
respectives (Cp^n) de chacun des sources sismiques (Si) aux frequences fondamentales et a 
leurs harmoniques, apr&s completion d'un cycle de mesures, les sismogrammes que Ton 
peut associer aux operations de surveillance sismique active ; et 

5 f) on forme par inversion dans le domaine des temps, k partir de la partie (P Ptn ) 

relevant de Tactivite microsismique passive, les signaux microsismiques sous-jacents 
contenus dans les enregistrements. 

8) Methode selon la revendication 7, caract£ris6e en ce que la contribution spectrale 
(C P xu) est obtenue en multipliant une fonction de transfert (T p ,j, nir ) entre une ondelette 

10 caract6ristique de la source et un sismogramme associe au r£cepteur r, par une ondelette 
caract6ristique de la source. 

9) M6thode selon la revendication 8, caracteris£e en ce que Ton met 
continuellement k jour la dite fonction de transfert. 

10) M6thode selon la revendication 9, caracterisee en ce que la mise k jour de la dite 
15 fonction de transfert (T Pt i tI1 , r ) est obtenue au cours d ! un cycle courant k partir d f une 

estimation (T Pt i tI1 _ 3i r)faite au cours d'un cycle pr£c6dent et d'une estimation initiale (To p ,i, n , r 
) faite au cours du cycle courant par la relation : 



Tpj, n .r =(l-h)T p>u _ l>r +hTo pa>n>r 
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